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ABSTRACT 

We present spatially-resolved imaging obtained with the Australia Telescope Compact Array 
(ATCA) of three CO lines in two high-redshift gravitationally lensed dusty star-forming galaxies, 
discovered by the South Pole Telescope. Strong lensing allows us to probe the structure and dynamics 
of the molecular gas in these two objects, at 2 : = 2.78 and z = 5.66, with effective source-plane res¬ 
olution of less than 1 kpc. We model the lensed emission from multiple CO transitions and the dust 
continuum in a consistent manner, finding that the cold molecular gas as traced by low-J CO always 
has a larger half-light radius than the 870 ptm dust continuum emission. This size difference leads to 
up to 50% differences in the magnification factor for the cold gas compared to dust. In the z=2.78 
galaxy, these CO observations confirm that the background source is undergoing a major merger, 
while the velocity field of the other source is more complex. We use the ATCA CO observations 
and comparable resolution Atacama Large Millimeter/submillimeter Array dust continuum imaging 
of the same objects to constrain the CO-H 2 conversion factor with three different procedures, finding 
good agreement between the methods and values consistent with those found for rapidly star-forming 
systems. We discuss these galaxies in the context of the star formation - gas mass surface density 
relation, noting that the change in emitting area with observed CO transition must be accounted for 
when comparing high-redshift galaxies to their lower redshift counterparts. 

Subject headings: galaxies: high-redshift — galaxies: ISM — galaxies: star formation — ISM: 
molecules 
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1. INTRODUCTION 

Carbon monoxide hereafter CO) has long 

been known as a tracer of molecular hydrogen gas in 
galaxies. Molecular gas is the fuel for new generat ions 

and 
3un¬ 


of stars (for recent re views, see Bolatto et al. 2013 
Carilli fc Walter|2013 ), so accurately diagnosing its a 


dance, kinematics, and morphology can shed light on 
the astrophysics of star formation. The most intense 
bouts of star formation in the universe appear to oc¬ 
cur in dus ty, star-forming ga laxies (DSFGs) at high red- 
shift (e.g., Casey et al.|2014 ). These galaxies are heavily 
enshrouded in dust, which absorbs the ultraviolet radi¬ 
ation from massive young stars and reradiates at far- 
IR and submillimeter wavelengths. These galaxies lie in 
contrast to the bulk of the high-redshift galaxy popula- 


ally isolated systems (e.g., 

Forster Schreiber et al.| 

2009 

Daddi et al. 

2010 

Tacconi 

et al. 

2013). Such “normal' 


emission, which generally excludes the highly obscured 
DSFG population. Together with rapid star formation. 


gas (> 10^° Mm; e.g., 

Greve et al. 2005 

Ivison et al. 2011 

Bothwell et al.|2013a 

that make up a signihcant traction 


2010 Ivison et al.||2011 Carilli fc Walter||20ta ). 


The most extreme L)ShGs ar e likely to be galaxies 


undergoing major mergers (e.g., |Narayanan et al. 

20101 

IHayward et al. |2012| |Fu et al.| 

20131 llvison et al. 


though see |Garilli et al.||20101 

Hodge et al.||2012 

2015 


(SFRs) enhanced by gas being tunneled to the center 
of the system after being disrupted during the collision. 
The merger kinematically manifests as a disordered ve¬ 
locity field or multiple comp onents closely separated in 
posit ion and/or velocity (e.g., Engel et al.|2010 Fu et al. 


2013 1. Such an extreme level of star formati on can likely 
only be sustained for a period of < 100 Myr (Greve et al. 


20051, and thus the brightest DSFGs are also relatively 


rare. 


Extensive effort has gone towards studying gas and 
dust both in the local universe and at high redshift in or¬ 
der to understand the physics and history of star forma¬ 
tion. One of the most studied correlations is the power- 
law relationship between the gas s urface density, Sp-as, 
and the SFR su rface density, Ssfr ( Schmidt||1959 Ken- 
|nicutt|1998[ see |Kennicutt fc Evans||2012| tor a recent re- 
view). On scales of a few hundred p arsecs, the two quan¬ 


tities appear line arly related (e.g., Schruba et al. 2011. 
Leroy et al.||2013), though power-law exponents ranging 


from sub-linear to quadratic have also been theoretically 
predicted and obs ervationally confirmed depending on 
methodology (e.K.,|Krumholz et al.|2009l |Liu et al. 12011 
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Faucher-Giguere et al. 2013"Shetty et al. 2013). I'his 
star formation (SF) relation (or Schmidt-Kennicutt rela¬ 
tion) is one ingredient in many theoretical prescriptions 
for star formation, so understanding its mathematical 
form and range of applicability is important for under¬ 
standing the buildup of stellar mass. 

The steps to derive a molecular gas mass from the lu¬ 
minosity of a low-J GO transition are not straightfor¬ 
ward, and a variety of techniques have bee n presented 
in the literature (e.g., Bolatto et al. 2013). The con¬ 


version factor, acOi varies with the kinematic state of 
the gas (through the escape fraction of GO photons) and 
the gas metallicity (through CO formation and destruc¬ 
tion processes). A variety of observations suggest that 
a value of aco 3.6 — 4.5 Mq pc“^ (Kkms“^)“^ (in¬ 
cluding a 36% mass contribution from the cosmological 
abundance of helium; hereafter we suppress the units of 
aco) is applicable to the Milky Way and nearby qui- 


metallicity (e.g., 

Solomon et al. 

1987 

Abdo et al. 

2010 

Sandstrom et al. 

2U13p. In regions o 

vigorous star tor- 


(e.g., Downes fc Solomon||1998 Tacconi et al. 2008). 


As a further complication, the high-EsFR galaxies that 
are much more common at high redshift are rare in the 
local universe, which makes their exploration more diffi¬ 
cult. Due to the faintness of the lowest transitions of 
CO, most high-redshift studies of molecular gas have 
used either spatially unresolved observations, or brighter, 
higher-J transitions with higher excitation conditions 
than the ground state, or in some cases both. An ad¬ 
ditional conversion from the observed CO transition to 
CO(l-O) is required, which depends on the temperature, 
density, and structure of the interstellar medium (ISM). 
Resolved observations of the lowest CO transitions are 
needed to test the SF relation on sub-galactic scales. 

Such high-resolution studies are aided by the use of 
gravitational leasing, in which a background object is 
magnified by a foreground structure, usually a massive 
elliptical ga laxy or galaxy cluster. For example, |Rawle| 
et al.| (|2014|) use high-resolution maps of dust continuum 
emission and [CII] and CO(1-0) emission to spatially and 
spectrally decompose a source at z = 5.2 predominantly 
lensed by a z = 0.63 galaxy, finding variations in the 
efficiency of star formation of a factor of ~6 x within a 


4kpc region in the source plane. Thomson et al.|(2015l 
probe ^ 100 pc scales in the z = 2.3 gal a xy SM M 


im- 


0102 (“the Eyelash”; Swinbank et al. ( 2010[ )). These 


authors studied the schmidt-Kennicutt relation in in¬ 
dividual star-forming clumps in this galaxy, and found 
evidence that the clumps are offset towards higher star 
formation efficiency compared to the galaxy as a whole. 


conducted 

by the South Pole Telsecope (SPT; 

Carlstrom 

et al. 

2011 

1 Vieira et al. 

2010 

Mocanu et al.| 

20131) and 

Herschel (J 

Negrello et al. 

2010 

Wardlow et al.|2013) have 


discovered large numbers of lensed DSFGs. Subsequent 
spectroscopy and high-resolution imaging have conhrmed 


redshifts and are lensed ( 

Weifi et al. 2013 

IVieira et al. 

20131 IHezaveh et al.||2013 

(Harris et al. 1201 

2 (Bussmann 

et al. 2013). In particular. 

Weill et al. 

(|2013 

1 used the At- 


to conduct a redshift survey of 26 DSFGs discovered by 
the SPT, finding a niedian redshift (z) = 3.5. Addition¬ 
ally, Aravena et al. (2013) and Aravena et ah, in prep. 
surveyed 18 ot these galaxies in low-J CO (either CO(l- 
0) or CO(2-l)). 

In this paper, we present high-resolution observations 
of low-J CO emission in two DSFGs from the SPT sam¬ 
ple performed with the Australia Telescope Compact Ar¬ 
ray (ATCA). Both objects have been observed at compa¬ 
rable (^0.5”) resolution by ALMA at 870/rm, with lens 
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models determined from these data (Hezaveh et al. 2013). 
SPT-S J053816-5030.8 (SPT0538-50), at z = 2.78, is rep- 
resentative of the typical DSFG population in redshift, 
870 ^m flux density, and dust temperatur e . This object 
was studied in detail by |Bothwell et al.| (|2013b|), who 
showed evidence for two velocity components in CO(7- 
6) separated by Av ^ 350kms~^. IntriRu ingly, the lens 
model of this source (Hezaveh et al.||2013) also required 
two dust components to reproduce the ALMA data, sug¬ 
gesting a possible physical connection between the ve¬ 
locity structure and the continuum structure. SPT-S 
J034640-5204.9 (SPT0346-52), at z = 5.66, is among the 
highest-redshift DSFGs known. The ALMA lens model 
indicates that it is also the most intrinsically luminous 
object in the SPT sample, and its Ssfr approaches or 
surpasses th e Eddington limit for radiation pressure on 
dust grains (|Thompson et al.|2005|) . While not represen¬ 
tative of the typical USh'G in the SPT sample, SPT0346- 
52 allows us to study the conditions of the ISM at their 
most extreme. 

The layout of this paper is as follows. In Section we 
describe the ATCA 7 mm and 3 mm observations m CO 
lines in these two objects. In Section]^ we describe our 
procedure for modeling the effects of gravitational lensing 
in both the ATCA and ALMA data. The morphological 
and kinematic results of these lens models are given in 
Section]^ In Sectionwe use the lens modeling results 
to discuss the effects m preferential source magnification, 
determine the aco factor in each source, and place these 
sources in the context of the SF relation. We conclude 
in Section We adopt the WMAP9 ACDM cosmology, 
with (ilm, Ha, Hn) = (0.286, 0.713, 69.3kms“^ Mpc” 


(Hinshaw et al. 2013). Throughout, we define the total 
infrared luminosity, Lir, to be integr a ted ov er rest-frame 
8-1000 /im, and assume a Chabrier (2003) initial mass 
function. 


was observed for flux calibration; when this source was 
not available, the flux level of the bandpass and gain cal¬ 
ibrators from adjoining observing dates was used to set 
the amplitude scale. Repeated observations of amplitude 
calibration sources indicate that the absolute flux scale 
is accurate to within 10%. _ 

The two S PT DSFGs were also observed by |Aravena| 
et al. (2013|) and Aravena et ah, in prep, using the com- 
pact ATCA configurations H75 (SPT0538-50, project 
C2655) and H214 (SPT0346-52, project C2744) in 2012 
July and October, respectively. For our present purposes, 
the baselines provided by these compact array configu¬ 
rations provide sensitivity to extended emission and an 
estimate of the total flux. The integrated spectrum of 
each source derived from these data are shown in Figs.j^ 
and The spectrum o f SPT 0538-50 has two peaks, as 
noted by I Aravena et al. (2013), which we discuss further 
in Sectionl3 

The three array configurations provide good uv cover¬ 
age on baselines from ^100 m - 6 km, and were combined 
and inverted using natural weighting. We ensure proper 
normalization of the noise levels of each dataset by differ¬ 
encing successive pairs of visibilities on the same baseline 
and polarization. Naturally-weighted channel maps of 
each galaxy are shown in Fig. IT and the upper panels of 
Fig.[^ We show channel maps m CO(2-l) in SPT0346-52 
in 200kms~^-wide channels, and separate the CO(l-O) 
line of SPT0538-50 into the red and blue velocity compo¬ 
nents seen in the integrated spectrum (two channels, ap¬ 
proximately 350 kms“^ wide). Note that the red velocity 
component of SPT0538-50 is not significantly detected, 
as the weak line flux is spread over several synthesized 
beams. The Icr sensitivities of these maps are 54 /rJy 
beam“^ (170 mK) per 200kms“^ channel in SPT0346- 


52 and 36 /xJy beam ^ 
in SPT0538-50. 


(53 mK) per 350 km s ^ channel 


2. OBSERVATIONS 

Table □ summarizes the ATCA observations of CO 
lines in SPT0346-52 and SPT0538-50. Further details of 
these observations are given below. We also make use of 
ALMA 870 ^m cont inuum observat i ons of these objects, 
described further in I Hezaveh et al.l (12013 1. 


2.1. ATCA Observations: 7mm Band 

SPT0538-50 and SPT0346-52 were observed with the 
ATCA 7mm receivers in CO(l-O) and CO(2-l) at ob¬ 
served frequencies of 30.45 and 34.64 GHz, respectively, 
as part of project IDs C2892 and C2983. The sources 
were observed using a 6-km extended array configuration 
over the course of 12 nights in 2013 October-November 
and a compact 1.5-km array configuration over 6 nights 
in 2014 January and May. The two CO lines are red- 
shifted to similar frequencies and can be observed with¬ 
out retuning the two 1 GHz-wide basebands available us¬ 
ing the Compact Array Broadband Backend (CABB), 
allowing bandpass and absolute flux calibration to be 
shared between the two sources observed in a single 
track. The bright quasars PKS1921-293 and PKS0537- 
441 were observed for bandpass calibration, while the 
quasars PKS0322-403 and PKS0537-441 served as com¬ 
plex gain calibrators for SPT0346-52 and SPT0538-50, 
respectively. For most tracks, the quasar PKS1934-638 


2.2. ATCA Observations: 3mm Band 

For SPT0538-50, we also observed the CO(3-2) line, 
redshifted to 91.35 GHz, using the hybrid H168 array 
configuration on 2013 August 20 in project C2816. The 
other 1 GHz-wide baseband was tuned to 94 GHz. The 
quasar PKS0537-441 was again used for bandpass and 
complex gain calibration, while Uranus was observed 
for flux calibration. The absolute flux scale at 3 mm 
is expected to be accurate to within ~15%, again in¬ 
ferred from repeated observations of amplitude calibra¬ 
tion sources. 

The data were continuum-subtracted and imaged us¬ 
ing natural weighting to maximize sensitivity to weak 
emission, giving a synthesized beam of ^3.lx 2.2”. This 
resolution is sufficient to marginally resolve the source, as 
seen in the bottom panels of Fig.[^ where we have again 
imaged the line in each of the two CO velocity peaks 
separately. These maps reach a sensitivity of 420 fiJy 
beam“^ (8.6 mK) in each 350kms~^ channel. Due to 
the larger synthesized beam size and higher line flux of 
the CO(3-2) line compared to the CO(l-O) line, we are 
also able to detect the weak, red velocity component in 
these data, at ~ 6a significance. 

We also significantly detect the dust continuum emis¬ 
sion at 3.3 mm using the line-free channels of both base¬ 
bands, reaching a sensitivity of 90 /iJy beam“^. The dust 
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Figure 1. Spectrum and channel maps of CO(2-l) observed in SPT0346-52. Left: Integrated spectrum derived from short, unresolved 
baselines, as presented in Aravena et al., in prep. The four 200kms“^ channels we model are colored. The dashed line shows the intrinsic 
spectrum using the magnification factors in Table|^ multiplied by 4x for clarity. Right panels: Channel maps of the high-resolution ATCA 
CO(2-l) observations pre sented here, colore d as m the left panel. The greyscale image of the lens galaxy is from co-added HSTfWFCS 
F140W + F160W images | |Vieira et al.|2013t . Thin black contours show the ALMA 870 /rm dust continuum images. The ATCA images are 
shown in steps of 2cr starting at ±3f7 (Icr = 10.6 mjy kms~^ bearn”^). The ALMA images are shown in steps of 5cr (Icr = 0.80 mjy bearn”^). 
Both datasets reach a resolution of approximately 0.45” x 0.65”. 


continuum emission closely resembles the emission from 
the red portion o f th e CO(3-2) line, which we discuss 
further in Section lOl 


3. LENS MODELING 

To derive the intrinsic gas and dust properties of the 
two DSFGs presented here, we must quantify the effects 
of gravitational leasing. Our lens mod e ling p rocedure 
follows that described by Hezaveh et al. (2013). Briefly, 
the lens mass profile is represented by a Singular Isother¬ 
mal Ellipsoid (SIE). For SPT0346-52, the model also 
strongly favors the existence of an external shear com¬ 
ponent whose axis is aligned with another galaxy ~3” 
east of the primary lens. For both sources, the ALMA 
870 ^m data are of much higher significance than the 
ATCA observations presented here, so we use the best- 
fit lens properties derived from the ALMA continuum 
data to model the source-plane in the ATCA CO data. 

The lensed CO source is represented by a parame¬ 
terized model consisting of one (SPT0346-52) or two 
(SPT0538-50; see below) symmetric Gaussian light pro¬ 
files in each modeled velocity channel. Each profile has 
up to four free parameters, namely, the two-dimensional 
centroid of the source and its intrinsic flux and size. 
While this source-plane model is undoubtedly overly sim¬ 


plistic, it allows the derived properties of each velocity 
bin to be compared in a straightforward manner. Us¬ 
ing a parametric model additionally avoids over-fitting 
the data using a large number of free parameters, as 
in most pixel-based reconstruction techniques, which are 
more appropriate for very high resolution, very high S/N 
observations. 

ATGA and ALMA both measure the Fourier compo¬ 
nents (visibilities) of the sky at the two-dimensional spa¬ 
tial frequencies defined by pairs of antennas. Rather 
than comparing to reconstructed images, where there are 
strong correlations between pixels, we fit lensing models 
directl y to the measured visibilities. As in |Hezaveh et al.| 
(2013|), we use a Bayesian Markov Chain Monte Garlo 
(iVlGMC) fitting procedure. At each MCMC step, we 
generate a model lensed image from a given set of source 
parameters. We then invert this image to the Fourier do¬ 
main and interpolate the model visibilities to the mea¬ 
sured uv coordinates of the ATCA data, using the 
metric to determine the quality of the fit. 

Lens modeling of the 870 /rm dust conti nuum emission 
of both sources was previously presented in |Hezaveh et al.| 
(2013) using ALMA data with approximately 1.5” reso- 
lutioii. In the present work, we additionally make use of 
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Figure 2. Spectra and channel maps of CO(l-O) and CO(3-2) observed by ATCA in SPT0538-50. Left panels: Integrated spectra of each 
line deriv ed from short baselines i n which the source is unresolved. We separate the two velocity components seen in these and other spectral 
lines as in|Bothwell et al.|l|2013b||. Top center and right panels: Channel maps of CO(l—0) of each of the colored velocity bins in the top left 
panel. The greyscale image is a co-added HST /WFC3 F140W + F160W image. Also shown as thin black contours are the ALMA 870 pm 
dust continuum observations. ATCA contours are shown in steps of 2 (t, starting at ±3(7 (Icr = 12.7mjy kms“^ beam“^), with a resolution 
of approximately 0.7” X 1.2”. The red velocity component is not significantly detected at the depth and resolution of these observations. 
ALMA contours are shown in steps of lOcr starting at +5 (t (Itr ~ 0.46 mjybeam“^). Bottom center and right panels: Channel maps of 
CO(3-2) in each of the colored channels in the lower left panel. Thick black contours show the 3.3mm continuum emission from the same 
data. All ATCA contours in these panels are shown in steps of 2a starting at ±3(7. For the CO(3-2) data, la = 145 mjykms“^ beam“^; 
for the 3.3mm continuum data, la ~ 70 /rjybeam“^. The ATCA data have a resolution of approximately 2.2”x3.3”. It is clear from this 
figure that the red CO(3-2) component is morphologically similar to the dust emission at the same wavelength, dominated by the bright, 
compact dust component found in the ALMA lens models (see text). 


higher resolution data taken as part of the same ALMA 
program, but which were not yet available at the time of 
publication of Hezaveh et ah. These new data were taken 
in an extended array configuration available in ALMA 
Cycle 0, and reach ~0.5” resolution. Further o bserva¬ 
tional details of this ALMA program are given in |Heza-| 


veh et al. (2013), while the updated 870^m lens models 
will be presented in full in Spilker et ah, in prep. For 
both sources, the updated lens models using the higher- 
resolution data are qualitatively and quantitatively sim¬ 
ilar to those derived by Hezaveh et al. Source properties 
relevant to this work are summarized in Table [21 

In SPT0346-52, the background source is well-fit by 
a single elliptical Sersic light profile with a circularized 
half-light radius of 610 ± 30 pc. The source’s intrinsic 
870 ^m flux density of 19.6 ± 0.5 mJy is magnified by a 
factor of ft = 5-5 ± 0.1. 

In SPT0538-50, the updated lens model requires two 
source-plane components to fit the data, as in Heza¬ 
veh et al. The source consists of a faint, diffuse dust 


component of intrinsic flux density 1.5 ± 0.2 mJy, half- 
light radius 1.25 ± 0.14 kpc magnified by ^ = 23.1 ± 2.2, 
and a brighter, more compact component of flux density 
3.7 ± 0.4 mJy and half-light radius 470 ± 50 pc magnified 
hy ft = 18.9 ±2.2. 


4. RESULTS 
4.1. SPT0346-52 

As shown in Fig. we detect the CO(2-1) emission 
from SPT0346-52 at peak significance of >5cr in four 
consecutive 200kms~^ wide channels. A single sym¬ 
metric Gaussian source-plane component for each fre¬ 
quency channel is sufficient to model the observed emis¬ 
sion with residuals consistent with noise in all cases. 
The derived source-plane structure is shown in Fig. 
where we have truncated the models of each channel at 
the half-light radius for clarity. The line emission blue- 
ward of ~ — lOOkms”^ is significantly offset from the 
redder emission. We have modeled alternative channel¬ 
izations of the data, with velocity bins ranging from 100- 
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Table 1 

Observational Summary 


Source Line 

^obs (GHz) 

Dates 

Array Configuration 

Time On-Source (h) 

Beam Size^ 

<7^ (/xjy) 

SPT0346-52 CO(2-l) 

34.636 

2012 10/04-10/09 

H214 

5.5 

6.4” X4.8” 

142 



2014 01/17-01/18 

1.5B 

6.0 

1.0”x0.6” 

180 



2014 05/10-05/17 

1.5D 

7.2 

1.6”x0.8” 

178 



2013 10/23-11/08 

6A 

26.1 

0.5”x0.3” 

66 

SPT0538-50 CO(l-O) 

30.450 

2012 07/28-07/29 

H75 

8.5 

19”X13” 

163 



2014 01/18-01/20 

1.5B 

11.9 

0.9”x0.5” 

119 



2014 05/15-05/17 

1.5D 

11.3 

1.2”x0.6” 

96 



2013 10/24-11/08 

6A 

45.7 

0.5”x0.3” 

42 

SPT0538-50 CO(3-2) 

91.345 

2013 08/20 

H168 

8.7 

3.1”X2.2” 

420 

^ Beam size for a naturally- 

weighted image. 





° rms sensitivity in 200 km s ^ (SPT0346-52) or 350 km s ^ (SPT0538-50) channels. 





Table 2 

CO and Dust Lens Modeling Results 

Source 

Line 

Component 

CO Lens Model Properties 

MCO Leo 

(lO^'^ Kkms“^ pc^) 

reff,CO 

(kpc) 

log(Egas)"*^ 

(Mq pc 2) 

SPT0346-52 

CO(2-l) 

-400kms“^ 

5.0 ± 0.6 

1.6 ± 0.2 

0.70 ± 0.13 

3.82 ± 0.32 



-200kms“^ 

5.8 ± 0.7 

2.3 ± 0.3 

1.86 ± 0.28 

3.14 ± 0.27 



0kms“^ 

10.0 ± 0.1 

1.5 ± 0.1 

0.77 ± 0.07 

3.72 ± 0.26 



+200kms“l 

10.1 ± 0.3 

1.4 ± 0.1 

0.73 ± 0.09 

3.74 ± 0.27 

SPT0538-50B 

CO(l-O) 

blue 

15.7 ± 2.3 

1.7 ± 0.3 

2.34 ± 0.43 

2.86 ± 0.26 

SPT0538-50A 


red 

- 

- 

- 

- 

SPT0538-50B 

CO(3-2) 

blue 

25.7 ± 1.2 

0.94 ± 0.13 

1.03 ± 0.14 

- 

SPT0538-50A 


red 

22.6 ± 1.2 

0.51 ± 0.07 

0.68 ± 0.12 

- 





Dust Lens Model Properties 



Source 

Component 

M870 Aim 

SsTO //m 

Lia 

^eff,870 fj,m 

log(SsFR)l> 




(mjy) 

1012 Lo 

(kpc) 

(M© yr-lkpc-2) 

SPT0346-52 


5.5 ± 0.1 

19.6 ± 0.5 

36.3 ± 5.4 

0.61 ± 0.03 

3.19 ± 0.08 

SPT0538-50A 

bright 

18.9 ± 2.2 

3.7 ± 0.4 

3.51 ± 0.50 

0.47 ± 0.05 

2.40 ± 0.11 

SPT0538-50 B 

faint 

23.1 ± 2.2 

1.5 ± 0.2 

1.48 ± 0.24 

1.25 ± 0.14 

1.18 ± 0.12 


Note. — When modeling the CO emission of both objects, we fix the parameters of the lens to their best-fit values derived from the 
ALMA data. For SPT0538-50, we also fix the location of the A and B source components to the positions derived from the ALMA data. 
^ Assuming the average CO-H 2 conversion factors determined in Section |5.2[ 1.4 for SPT0346-52 and 1.5 for SPT0538-50. We propagate 
an uncertainty of 50% into the determinations of Sgas- 
^ Surface densities determined within reff. 


400kms“^, and all channel widths point towards the 
same overall structure. Derived source properties are 
given in Table 

The structure seen in Fig. is difficult to interpret. At 
the depth of these data, the -mocity structure is not obvi¬ 
ously consistent but not clearly inconsistent with large- 
scale disc rotation or other bulk motion. As we have 
only modeled four consecutive velocity channels, it is dif¬ 
ficult to rule out either ordered or disordered kinematics. 
Taking the best-fit centroid of each velocity component 
at face value, the molecular gas in SPT0346-52 appears 
more consistent with a merging system than with the 
massive rotating discs seen in normal star-forming galax- 
ies at mod erate redshift by, e.g., |F5rster Schreiber et al.| 
(2009) and Tacconi et al. (2013)' If, on the other hand, 


the centroid of the -200 kms~^ component is, in fact, be¬ 
tween the bluer and redder channels, a position-velocity 
curve resembling a rotating disc could result. Similar 


positio n-ve locity diagrams wer e found by |Riechers et al.j 
(2008) and Deane et al. (2013) in the source-plane struc- 
ture of lensed quasars at z = 4.1 and z = 2.3, using 
data of similar significance to that presented here. These 
authors interpreted their data as suggestive of rotation. 
Deeper observations are necessary to draw stronger con¬ 
clusions about the velocity structure of SPT0346-52, as 
this would allow lens modeling of narrower velocity bins. 

The half-light radius of the background source is larger 
in low-J CO emission than in rest-frame 130 /rm contin¬ 
uum emission, and the 130 /rm emission appears closely 
associated with the blue half of the CO (2-1) emission. 
This tentatively suggests that the star formation in 
SPT0346-52 is proceeding in a compact region embed¬ 
ded within a larger reservoir of molecular gas. Alterna¬ 
tively, the CO(2-l) emission may trace a larger region of 
the galaxy owing to its higher optical depth. To test 
the degree to which a source with the same size and 
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Figure 3. Source-plane reconstruction of SPT0346-52. Each 
C0(2-l) velocity channel modeled is represented by a colored disk, 
truncated at the half-light radius. Colored “x” symbols and el¬ 
lipses show the centroid and la positional uncertainty on the cen¬ 
troid for each channel. The FWHM and location of the 870 ^m 
dust emission is shown as a thick dashed line, while the lensing 
caustic is shown as a thin black line. 


location as either of the CO-emitting components at 0 
and +200km/s could contribute to the 130 contin¬ 
uum emission, we re-fit the ALMA data with two source- 
plane components. We fix the size and position of each 
source, leaving only the flux density of each source as 
a free parameter. The size and position of one source 
is fixed to that derived from the ALMA data, while the 
other is fixed to the size and position of either of the 
two red CO channels. This is effectively a null test to 
determine how much flux density at rest 130 fim could 
be emitted from the same region as the 0 or -1-200kms~^ 
CO emission. This test indicates that a source co-located 
with either of the two reddest modeled CO channels con¬ 
tributes < 3% (Ict upper limit) of the total unlensed 
flux at 130 /im. This is also the fraction of Lir and, by 
proxy, SFR, that could arise from these locations under 
the assumptions of a uniform dust temperature across 
the source and no contribution from dust heating due to 
AGN activity. This limit is nevertheless consistent with 
CO/Lir ratios seen in local ULIRGs and z > 0.4 main- 
sequence galaxies (e.g., Ivison et al.||2011 Aravena et al. 
in prep.). 

We note that our finding of a large CO spatial extent 
compared to dust continuum is unlikely to be explained 
by the effects of interferometric filtering. While only a 
limited range of spatial frequencies are probed by both 
the ATCA and ALMA data, the data probe a similar 
range of radii in the uv plane. Additionally, as we fit 
directly to the visibilities, our model natively reproduces 
this filtering. The models presented here recover 85- 
105% of the flux observed in the most compact array 
configuration data. These data are included in the mod¬ 
eling, so it is unsurprising that the total flux should be 
recovered well. 


4.2. SPT0538-50 



-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 

Ax (arcsec) 


Figure 4. Source-plane reconstruction of SPT0538-50. As before, 
each component is truncated at its FWHM. In both panels, the 
bright/compact and faint/extended rest-frame 230/im dust com¬ 
ponents are represented by small and large dashed circles. The 
best-fit extent of the blue CO(l-O) velocity component is shown in 
the top panel, while the bottom shows the extent of both the blue 
and red CO(3—2) velocity components. Note that we are unable to 
model the red CO(l-O) velocity component, and that the centroids 
of each CO component were fixed to the best-fit positions derived 
from the ALMA continuum data. 


The integrated spectrum of this source shows two ve¬ 
locity components in CO(l-O) emission separated by 
350km/s (Fig. [Aravena et al.|j2013 Bothwell et al. 


2013b| ), and so we model the eniission from each ve¬ 


locity component separately. Unfortunately, our high- 
resolution CO (1-0) maps are too shallow to detect the 
faint red velocity component. However, this source was 
also observed in CO(3-2) emission, marginally resolving 
the source and significantly detecting both velocity com¬ 
ponents. These observations also detected the observed- 
frame 3.3 mm continuum emission at > 5cr. Here, we 
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discuss the results of the lensing inversion of these three 
datasets in turn. As before, we hx the parameters of the 
lens to those derived from the ALMA 870 /xm data. 

4.2.1. Dust Continuum Emission 

As p reviously described and as found by |Hezaveh et al.| 
(j2013l), we require two source-plane components to ht the 
870 ^m continuum emission observed by ALMA (rest- 
frame 230/im). The source consists of a faint, diffuse 
dust component of intrinsic flux density 1.4±0.4mJy, 
half-light radius 1.2 ± 0.3 kpc magnified by /x = 24 ± 4.5, 
and a brighter, more compact component of flux density 
3.5 ±0.7 mJy and half-light radius 460 ± 90 pc magnified 
by /i = 20 ± 4. 

We have now also detected dust continuum emission at 
observed-frame 3.3 mm (rest-frame 870 /xm) in the line- 
free channels of our ATCA 3 mm data. As the continuum 
is much more weakly detected in the ATCA data than 
the ALMA data (Fig. [^, and as dust gives rise to the 
observed emission at both wavelengths, we fit the ATCA 
data with two source-plane components with positions 
and sizes fixed to the best-fit values derived from the 
ALMA data. This leaves only the fluxes of each compo¬ 
nent as free parameters. 

Assuming the source-plane morphology of the dust 
emission is the same at both rest-frame 230 and 870 /xm, 
we derive intrinsic fluxes of 55 ± 15 /xJy for the bright, 
compact (A) component and 5 ± 9 /xJy for the faint, 
diffuse (B) component. By calculating the flux ratio 
of the components at each MCMC step, we rule out 
the possibility that the two components have the same 
2.5:1 flux ratio seen at 870/xm at the ^ 2.2cr level, 
tentative evidence that the two components have dif¬ 
ferent metallicities, dust temperatures and/or opacities. 
Further multi-wavelength high-resolution continuum ob¬ 
servations could resolve this issue, though the atmo¬ 
sphere limits such prospects at wavelengths shorter than 
observed-frame 350 /tm (rest-frame 92 /tm). 

4.2.2. CO(3-2) Emission 

In star-forming galaxies, the CO(3-2) transition is sig¬ 
nificantly more luminous than lower-lying transitions; for 
thermalized level populations, the integrated line flux 
scales as J^. This, combined with the larger beam size 
of our ATCA 3 mm observations, allowed us to detect 
the fainter red velocity component in CO(3-2) emission 
while this component remained undetected in CO(l-O). 

As clearly seen in the images of Fig. the spatial dis¬ 
tributions of the blue and red CO(3-2)emission appear 
significantly different. The red velocity component ap¬ 
pears similar to the rest-frame 870 /xm continuum emis¬ 
sion simultaneously observed by ATCA. This emission is 
dominated by the bright, compact (A) dust component, 
as described in the previous section. We begin exploring 
the lensing inversion by htting a single Gaussian profile 
to each velocity component, allowing the position, size, 
and flux of each source to vary. The best-ht position of 
the blue velocity component is within the Itr uncertain¬ 
ties of the location of the diffuse B dust component seen 
in the ALMA data, and inconsistent with the location 
of the compact A component. While the positional un¬ 
certainties are large (~0.13”), the best-ht position of the 
red velocity component is consistent with the location of 
the bright, compact dust component seen in the ALMA 


maps. This suggests that the two velocity components 
are in fact associated with the two dust components. 

Motivated by this association, we re-£t the CO(3-2) 
data of both the red and blue velocity components with 
two source-plane components each, hxing the positions 
of each component to the best-ht positions of the dust 
components in the continuum lens model, allowing the 
source size to remain a free parameter. This is a test 
to determine whether the CO emission of each velocity 
component can be uniquely associated with one of the 
two dust components. The results of this test strengthen 
the hypothesis that each velocity component is associ¬ 
ated with only one of the two dust components - for 
each velocity component, the lens modeling prefers that 
only one of the two source-plane components have pos¬ 
itive hux and non-inhnite size. Again, we hnd that the 
brighter blue velocity component is associated with the 
faint, diffuse (B) dust component, and the fainter red ve¬ 
locity component is associated with the bright, compact 
(A) dust component. 

Figure shows the results of modeling the CO(3-2) 
emission, where we have hxed the locations of the blue 
and red velocity components to the locations of the faint 
and bright dust components, respectively. For each ve¬ 
locity channel, we allow the hux and size of the modeled 
source-plane component to vary. The models imply that 
the CO(3-2)-emitting molecular gas has approximately 
the same extent as the dust emission in both components. 

Having associated the two dust components seen in 
the lensing reconstruction with the two separate veloc¬ 
ity components seen in the integrated CO line spectra, 
we strengthen the argument that SPT0538-50 is indeed 
a pair of merg ing galaxies, as also posited by Bothwellj 


et al. (2013b). Their arguments, based on the high 
yi'K surface density, high specihc SFR, and suppressed 
hne structure lines of SPT0538-50, are conhrmed based 
on our high-resolution kinematic observations. The two 
merging galaxies are separated by 1.3 kpc in projection 
and ^ 375kms~^ in velocity. This is comparable to 
what is seen in the local ULIRG Arp220, a late-stage 
merger which shows two nuclei sep arated by approxi¬ 
mately 400 pc and ^ 2 50 km s~^ (e.g. 


Sakamoto et al.|[2009|. In contrast, S 


ppr 

Scoville et al. 1997 
^’T’U538-50 appears 


to be in a more compact merger than the z = 5.24 Her- 
schel-selecte d lensed galaxy HLS0918 studied by |Rawl^ 


et al. (2014), which consists of four spectral components 
separated by 4 kpc and ~ 840kms“^. SPT0538-50 ap¬ 
pears to conform with the idea that most DSFGs reach 
their extreme SFRs through major merger activity (e.g., 
Engel et al.|[^ |. 


4.3. CO(1-0) Emission 

As seen in Fig. our observations of CO(l-O) are in¬ 
sufficiently deep to detect the faint, red velocity compo¬ 
nent in our high-resolution imaging data. We do, how¬ 
ever, clearly detect the brighter blue velocity component 
at comparable resolution to the ALMA 870 /xm data. The 
morphology of the CO(l-O) emission is clearly different 
from that of the dust continuum emission, which is dom¬ 
inated by the bright, compact dust component described 
previously. While we cannot model the red hne compo¬ 
nent in CO(l-O), we proceed by discussing the blue line 
component, which can be modeled. 
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Having established in the previous section that the blue 
line component is spatially associated with the faint, dif¬ 
fuse dust component, we model the blue CO(l-O) emis¬ 
sion with a circularly-symmetric Gaussian source-plane 
component with position fixed to that derived from the 
ALMA continuum data. The free parameters are the 
source flux and size, as in our models of the CO(3-2 
emission. The best-fit model source, shown in Fig 
and described in Table leaves residuals consistent wit 
noise. Comparison of the CO(3-2) and CO(l-O) intrin¬ 
sic line fluxes implies a CO brightness temperature ratio 
of r^i ~ 0.6, similar to the ratios determined for other 


tne 

i 

ith 


age ratio for the SPT DSFGs themselves (e.g., Danielson 

et al. 

*M' ' 1 ' 

2011| Bothwell et al. 

2013a 

Spifker et ai. 

2014). 


find that the CO (1-0) emission in SPT0538-50 is signif¬ 
icantly extended compared to the dust emission. Again, 
this may indicate that the star formation in this system 
is proceeding in a compact region embedded in a larger 
reservoir of molecular gas. We return to this discussion 
in the next section. 


5. DISCUSSION 

5.1. Source Sizes and Differential Magnification 

As can be seen in Figures and we find that the 
molecular gas reservoirs traced by low-J CO transitions 
have larger half-light radii than the emission from the 
dust continuum. If we attribute all of the dust emis¬ 
sion to star formation (as opposed to, for example, dust 
heated by AGN activity), this implies that the intense 
star formation in these galaxies is limited to relatively 
small regions embedded in much larger reservoirs of 
molecular gas. The kinematics of both galaxies are plau¬ 
sibly consistent with the disruption of secular rotation, 
causing vast amounts of molecular gas to fall towards 
dense, compact star-forming regions. 

In SPT0538-50, we find that the CO(l-O) is similarly 
extended compared to CO (3-2), by more than a factor 
of 2x in the blue velocity com ponent. Size dif f erence s 
of this level were reported by ||Riechers et al.| (|2011|), 
comparing the and CO(l-O) sizes of the lensed DSFC 
SMM J09431-I-4700. Larger physical extents for low-J 
CO emission were also inferred by those authors and |Ivi-| 
[son et aL]|2011| by comparing the CO line widths be¬ 
tween low- and mid-J CO transitions. In contrast with 
the CO(I-O) emission, we find the CO(3-2) emission in 
SPT0538-50 to have roughly equal half-light radius as 
the star formation traced by the dust continuum. This 
can be taken as evidence that the CO(3-2) emission is 
more directly associated with ongoing star formation in 
this galaxy, in agre ement with studies bo th locally (e.g., 
Wilson et al. 2009|and at high-redshift (Bothwell et al. 


2(Jll)'T’acconi et al. 20131 which find an approximately 
linear relationship between the CO(3-2) luminosity and 
SFR. Indeed, high gas excitation conditions (T^in ~ 50 K, 
1000 cm“^) are needed to achieve brightness tem- 


riH. 


perature ratios near unity. A similarity i n size between 
CO(3- 2) and stellar light was also seen by Tacconi et al. 
(2013), who used rest-frame R-band RST images to de- 


termine the extent of the star formation in a large sample 
of z ^ 1 — 2 normal star-forming galaxies. 

While measurements are few, in Fig.[^we plot the half¬ 


light radii of star formation (traced by the rest-frame 
UV, dust continuum emission, or both) and molecular 
gas reservoirs for local and high-redshift sources. We re¬ 
strict this comparison to those galaxies with measured 
sizes in CO transitions with Jup < 3, since the effective 
source size changes significantly as a function of observed 
transition. Physically large molecular gas reservoirs are 
common, with an average star formation area filling fac¬ 
tor of ~ 55%, similar to the differ ence in size between 
gas and SF in local gala xies seen by|Bigiel & Blitz|(|2012 1 
and Zahid et al. (2014). This difference m size can be a 
potential source of bias when calculating surface-density 
quantities using sizes derived at different wavelengths - 
for example, using the source sizes determined from a SF 
tracer to calculate Sgas would lead to an over-estimation 
of the average surface density by a factor of 80%. 



reff (CO), kpc 


Figure 5. Relative sizes of star formation and molecular gas in 
galaxies. CO source sizes are confined to those objects for which 
low-J resolved measurements have been made. Star formation 
sizes are derived either from rest-frame UV measurements, dust 
continuum emission, or both. The solid line indic ates rofffSFj = 
rpfffC O ). The high-redshift sample is drawn fr o m [Younger et al.| 


(|2008ll; IDaddi et al .l (I 2 OIOI 1 ; [ivison et ah l ( 20111); IFu et al. h 2Ui:^;b: 
llodge et al. ] (pir^ ; f^ter el al.| ( |2aT7| l; ]t‘u et al.l ll^Ul^: jlviscm 
et al. l[?Ui;j[|. I'he loca l galaxies are drawn rrom|Ke^n et ah iWnorlr 


Ivisonl 


^ ^ |Kegan et al.|lp0ui|l 

and |Lj^oy et al. JpoM, where we have scaled the <JU sizes m Kegan 
et al. to matcti tae distances in Leroy et al.. 


This difference in source size and structure at differ¬ 
ent emitting wavelengths leads to differential magnifi¬ 
cation. Although gravitational leasing is achromatic, 
galaxies appear morphologically different at different 
wavelengths, leading to wavelengt h-dependent variations 
in the leasi ng magnification (e.g ., Blain et al.||1999 Ser¬ 


jeant 2012 Hezaveh et al. 2012). In both sources stud¬ 


ied here, the large total extent of the CO emission rela¬ 
tive to the dust continuum emission and the proximity 
of the sources to the leasing caustics lead to differences 
in magnification between dust and molecular gas of up 
to ~ 50%. In SPT0346-52, our data indicate that the 
gas emitting at positive systemic velocities is magnified 
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by approximately a factor of 2 x more than either the 
rest-frame 130 /rm continuum or the gas at bluer sys¬ 
temic velocities. Such a magnification gradient may in 
fact be a common feature of observing lensed galaxies 
in spectral lines in wh ich galaxies are expec ted t o have 
large physic al extents. Riechers et al. (2008) and Deane 
et al. (2013) both find similar levels of ditlerential magm- 
fication as a function of velocity in observations of low-J 
CO emission in the lensed quasars PSS J2322-I-1944 and 
IRAS F10214-I-4724, respectively. These results indicate 
that differential magnification can cause significant dis¬ 
tortion of the global SEDs of lensed galaxies, reaffirming 
the need to account for its effects when comparing lensed 
and unlensed sources. 

5.2. The CO-H 2 Conversion Factor 

With the wealth of high-resolution data on these two 
sources, we are presented with the opportunity to con¬ 
strain the conversion factor between CO luminosity and 
molecular gas mass, ocO) through multiple means. This 
allows us to measure the extent to which different meth¬ 
ods of determining aco are consistent within a single 
system. Below, we briefly describe each technique as ap¬ 
plied to SPT0346-52 and SPT0538-50. In each case, we 
assume the mass in the form of atomic hydrogen gas is 
negligible, as seems appropriate for dense, highly star¬ 
forming systems. We also assume the CO(2-l) line ob¬ 
served in SPT0346-52 is thermalized and optically thick. 


with L'r 


CO(i-o) ~ -^CO(2-i)- Given the vigorous star for¬ 
mation proceeding in this warm, dense system and that 
Tcmb(-z = 5.7) = 18K > Eup_co( 2 -i)/fcB> this assump¬ 
tion is justified. This ratio is also justified both observa- 
tionally (s ee, e.g., Fig. 45 ofiCasey et al. ||2014|) and the¬ 
oretically ( Narayanan &: Krumholz|^14 ). The results of 
this section are summarized in l'ig'[^ 

5.2.1. Gas-to-Dust Ratio 

Gas and dust are widely observed to be well-mixed in 
galaxies, with dust com prising approximately 1 % of the 
mass in the ISM (e.g., Sandstrom et al. 2013 Draine 
et al. 2014). Thus, if the mass in dust (Md) and the 
gas-to-dust ratio (^gdr) can be estimated, aco can be 
determined simply as Md(5GDR/Lgo- 

Extensive Spitzer and Herschel observations of the 
dust continuum emission in galaxies have rapidly ad¬ 
vanced our understanding of (5 gdr and the nature of 
the dust emission. The metallicity of the ISM affects 
^GDRj so that low-metallicity systems have more gas per 
unit dust mass. Meanwhile, most far-infrared dust SEDs 
can be accurately modeled using only a small number 
of free parameters. One common method is to fit the 
dust SED with a single-temperature modified blackbody 
function. However, such single-temperature models gen¬ 
erally un derestimate the dust mass in galaxies by a factor 
of ~ 2x (Dunne & Eales 2001 Dale et al. 2012). Dust 
at a single temperature cannot simultaneously lit both 
the long- and short-wavelength sides of the SED, which 
has contributions from dust heated to a range of equi- 
librium temperatur es. This effect has be en explored by 
Dunne et al. (2000) and Dale et al. (2012), and indicates 
that a more sophisticated, multi-component approach is 
needed. Such app roac hes have been develope d by, e.g.. 


exposed to a power-law distributi on of starlight intensi- 
ties, while the MAGPHYS code of da Cunha et al. (2008) 
decomposes the dust emission into difterent physlcally- 
motivated temperature regimes. The Draine & Li (2007) 
models have been used to constrain aco m samples of 
normal and intensely star-forming 2 ^ 0.5 — 4 Herschel- 
selected st y-forming galaxies ob served by |Magdis et al.| 
(2012) and Magnelli et al. (2012). 

Calculation of the dust mass requires knowledge of the 
dust mass absorption coefficient, defined here as Ksrotim, 
the value of the coefficient at rest-wavelength 870/xm, 
in units of m^kg“^. The dust mass scales as 
Most estimates of ksto for dust in nearby. Milky Way¬ 


like galaxies are close to either 0.045 (Li &: Draine 


2001 


Draine & Li 20071 

Scoville et al. 2014) or 0.075 (|Us- 

senkopt & Henning] J 

i994||Uunne et al.|200U| |James et al. 


terence in dust mass of ^70%. In this work, we adopt 


the dust mo dels implemented in MAGPHYS (da Cunha 
et al. 2008 ), which uses the dust opacity coefficient ot 


Draine & Li 


( 20071) and da Cunha et al. (2008). The the- 


oretical models of Draine &: Li] 2007) assume the dust is 


Dunne et al. (2000), KsTO/xm = 0.077. This value for the 
mass absorption coefficient is commonly used for other 
high-redshift rapidly star-forming systems. 

Einally, to calculate the total gas mass from the dust 
mass, we require knowledge of the gas-to-dust ratio, 
^GDRj which is known to vary with metallicity. As we 
have essentially no constraints on the metallicities of 
the se systems, we adopt th e average ^gdr determined 
by Sandstrom et al. (2013) in a large sample of local 
star-forming galaxies with approximately solar metallici¬ 
ties, ^gdr = 72 with ~0.2dex scatter. We note that the 
derived aco conversion factors are linearly proportional 
to the assumed i5gdr, and (5gdr its elf likely varies ap - 
proximately linearly with metallicity (Leroy et al.|2011 ). 

We determine the dust masses of SPi’034B-5 2 aii^d 
SPT0 5 38-50 by fi tting to the photometry given in |Weifi| 
et al.| (|2013|) and |Bothwell et al.| (|2013bl) , respectively. 
In both cases, we assume that ditterential magnification 
of the dust emission is insignificant, and for SPT0538-50, 
we assume the flux at each wavelength is divided between 
the two components in the same ratio as at rest-frame 
230/xm. This yields dust masses of 2.1 ± 0.3 x 10® Mq 
for SPT0346-52 and 9.0 ± 1.3 and 3.7 ± 0.6 x 10® Mq 
for the compact and diffuse components of SPT0538- 
50, where the uncertainties are statistical only. Eor 
SPT0346-52, this dust mass and the source size derived 
from the ALMA lens model imply that the dust reaches 
unit optical depth by rest-frame ~300 /xm, longer than 
the canonical wavelength of ^100/xm. Lower optical 
depths are possible by raising the dust temperature or 
effective source size. The prevalence of this effect in sam¬ 
ples of lensed DSFGs will be explored in more detail by 
Spilker et ah, in prep. 

Using our adopted ^gdr and dividing by the sum of 
the intrinsic CO luminosities of each channel yields a 
measurement of aco = 2.2 ± 0.6 for SPT0346-52. Eor 
SPT0538-50, we calculate aco for the blue velocity com¬ 
ponent only, yielding aco = 1-7 ± 0.4. Note that the 
uncertainties in these calculations account only for the 
statistical errors in the dust SED fitting and lens mod¬ 
eling procedure, and neglect systematic uncertainties in 
^870 /im and (jcDR, which are of order 100%. 
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5.2.2. Dynamical Constraints 

The CO-H 2 conversion factor can also be constrained 
using estimates of the total dynamical mass of galaxies 
and a process of elimination - the molecular gas mass is 
the remainder after all other contributions to the dynam¬ 
ical mass have been subtracted (e.g., stars, dark matter, 
HI gas, dust, etc.). Assuming that stars, molecular gas, 
and dark matter make up the vast majority of the total 
mass, then, we have 

Mdyn = M* -I- acoLco + Mdm, (1) 


where each of these quantities is measured within the 
same effective radius defined by the extent of the CO 
emission. 

Measuring each of these quantities is fraught with 
assumptions and systematic uncertainties. The stellar 
mass of SPT 0538-50 has been estimated by [BothweH] 


et al. (2013bl as 3.3 ± 1.5 x 10^°M©, while Ma et ah, 
in prep, only place an upper limit on the stellar mass of 
SPT0346-52, based on SED fitting to optical through far- 
infrared photometry. These estimates assume the stellar 
light is magnified by the same factor as the dust emis¬ 
sion, which may be inaccura te by up to a factor of ~50% 
(see Section 5.1 above and Calanog et al. 2014). The 
dark matter content within the region traced by CO is 
highly uncertain, but we adopt a contribution of 2 5% for 
consistency with the literature (jDaddi et al.|2010 1 , based 
on observations of z ^ 1.5 — 2 disk galaxies. 

A crude estimate of the dynamical mass can be derived 
from 

Mdyn = (2) 


with AV the FWHM line width, R the source effec¬ 
tive radius, and the gravitational constant G. Here the 
pre-factor 7 accounts for the detailed geometry of the 
source, and is of order unity. As the geometry of the two 
objects studied here is only somewhat constrained, we 
adopt 7 = 1 for si mplicity. Literature values of 7 range 


from y ~ 0.3 (e.g.. 

Neri et al. 112003 Tacconi et al. 12008 

Daddi et al. 2010 

) to 7 = 1.2 (e.g., liothwell et al.|2013a) 


In the following, we quote only statistical uncertainties 
on the derived values of acOi ignoring the much larger 
systematic uncertainties on the stellar mass, dark matter 
fraction, and true source geometry. 

For SPT0538-50, using the spatial separation of the 
two components of 1.35 ± 0.16 kpc and a line width of 
490kms“^ (from the full line profile), the above equation 
yields Mdyn ~ 7-5 ± 0.9 x 10^° M©. Combined with the 
stellar mass estimate given above and a 25% dark matter 
fraction, this yields aco = 1-4 ± 0.8. If the stellar mass 
is entirely concentrated in the blue velocity component, 
the implied gas fraction of the blue component is /gas = 
Mgas/(Mgas + Mstar) ~ 40%. As the red component must 
also contain some gas, 40% is a lower limit to the gas 
fraction of the entire system. 

We estimate the dynamical mass of SPT0346-52 using 
the separation in space and velocity of the emission at 
-400kms~^ and -|-200kms~^. This yields R = 1.8 ± 
0.2kpc and Mdyn = 1-5 ± 0.2 x 10^^ Mq. The stellar 
mass upper limit calculated by Ma et ah, in prep, is 
a factor of several higher than the dynamical mass we 
have estimated. To constrain ocoj we instead assume a 
range of gas fractions of /gas = 0.3 — 0.8. This range has 


been observed by 

Tacconi et al. 

(2013|) at lower redshifts 

(z ^ 1 — 2 ), and 

appears to evolve slowly with redshitt 

at z > 2 (e.g.. 

Bothwell et al. 

2013a Aravena et ah, 


in prep.). This range in gas fraction leads to a range 
of aco = 0-5 ~ 1-3; we adopt aco = 0-9 ± 0.5, where 
the uncertainty reflects only the range of /gas and the 
statistical uncertainty on Mdyn. 

Given the large systematic uncertainties inherent in 
each step of these calculations, we estimate that the un¬ 
certainty in these derived conversion factors is at least a 
factor of 2 X. Improving these estimates would require 
extensive high resolution multi-wavelength observations 
in order to better constrain the dynamical and stellar 
masses. Refined estimates of the dark matter contribu¬ 
tion would be even more challenging, and the best option 
may simply be to use the dark matter content derived 
from hydrodynamical simulations of galaxies. 


5.2.3. CO Luminosity Surface Density 

A third est imator of the CO-H^ co nversion factor was 
developed by|Narayanan et al.|(|2012D, based on hydrody¬ 
namical simulations coupled with dust and line radiative 
transfer. Those authors developed a suite of simulated 
galaxies in isolated and merging systems and provided 
a fitting formula to aco that depends solely on the CO 
line intensity and metallicity: 


aco = 10.7 X S]co3V(^/^0)°•®^ (3) 

where Sco is the CO luminosity surface density in units 
of Kkm/s. As in the dust-to-gas ratio method, we again 
assume solar metallicity for both DSFGs considered here; 
lowering the metallicity to half solar would increase these 
estimates by approximately 50%. For SPT0346-52, we 
average the values of aco determined for each modeled 
channel to find a value of aco = 0.78 ± 0.11. As we 
could not model the red velocity component in SPT0538- 
50, for this galaxy we calculate aco for the blue velocity 
component only, yielding aco = 1-4 ± 0 . 2 . 


5.2.4. Summary of CO-H 2 Conversion Factor 
Measurements 

We have constrained the CO-H 2 conversion factor 
in the two objects presented here using three indepen¬ 
dent methods. The three methods show reasonably 
good agreement with each other, with average values of 
aco = 1-3 for SPT0346-52 and aco = 1-5 for SPT0538- 
50. The uncertainty in these estimates is dominated by 
systematic, rather than statistical, errors of ^ 100 % for 
each method. Given the large systematic uncertainties 
inherent to each method, the general agreement between 
the three techniques is encouraging. Resolved CO and 
dust continuum observations of a larger sample of objects 
could reveal systematic differences between the various 
methods. 

For both objects, the derived conversion factors are 
similar to those determined for other rapidly star-forming 
objects, exhibiting low values of aco similar to most DS¬ 
FGs and unlike more quiescently star-forming objects. In 
Figure we place our measurements in the context of 
other studies that have also constrained the conversion 
factor in a wide variety of galaxies at z > 1. No clear bi¬ 
modality between “Milky Way-like” and “ULIRG-like” 
values is seen, as would be expected from a heteroge- 
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neous collection of galaxies with varying ISM proper¬ 
ties. Indeed, the sample used to d etermine the canon¬ 
ical “U LIRG-like” value of aco by [Downes &: Solomon 


(19981 also showed a fair a mount of variation. This lack 
ot birnodality (as noted by Narayanan et al.j 2012 and 
discussed in the next section] can influence the form of 
the Schmidt-Kennicutt star formation relation. 



Figure 6. Summary of constraints on the CO-H 2 conversion 
factor at z > 1 from the literature. Objects are color-coded by the 
method used to constrain aco shape-coded by type of object. 
SPT0346-52 and SPT0538-50 both appear to follow the general 
trend of decreasing aco with increasing Lir. Where applicable, 
Lir has been corrected for lensing magnification. If not given in 
the original references, measurement uncertainties for the literature 
objects are shown as 50%. Objects with constraints from multiple 
methods are shown multi ple times at the same Lit? . Literatur e 
objects are compiled fromIDaddi et a l. ||20 1 0||;|Ivison et al. | (|2011|>; 
Magdis et al.| 2011[l:[Swinbank et al.|(|2 )ll||;|ku et al. (|2Ui2[ | ::piodge 


et al.l ; Magdis et al. | [Magnelli et al.l 2U12 ; I Walter 

et al. ; lueane et al .l ||2Ui3b; fku et alTl (I 2 U I 3 V; IhLodge et al. 

'^^^^ne^afrlll^^Sl; iMessia^^^ TIM-' '- 




5 . 3 . The Star Formation Relation at z > 2.5 

We now turn our attention to the Schmidt-Kennicutt 
EsFR“Sgas relation. Figure presents samples of low- 
and high-redshift galaxies for comparison to the targets 
of this work. Our lens modeling and FIR SED provide 
measurements of the total star formation rate, molecular 
gas mass (inferred from low-J CO measurements and the 
aco measurements of the previous section), and sizes of 
the star-forming and molecular gas regions. The Esfr 
and Sgas values for these two objects show them to lie 
at the upper edge of the distribution of galaxies in this 
plane. 

To interpret the offset between the properties of the 
SPT galaxies and the larger sample consisting of both 
normal and more highly star-forming galaxies, we must 
consider the measurements that are required to con¬ 
struct such a plot. The high-redshift comparison galax¬ 
ies in Fig. [^are divided into two samples depending on 
whether their molecular gas sizes are determined from 


low-J CO (orange diamonds) or from other measure¬ 
ments, including star formation tracers and higher-J CO 
(purple squares). The galaxies with low-J CO size mea¬ 
surements, which provide determinations of Sgas that are 
most similar to those made at low redshift, are more 
consistent with the SPT sources than the other high- 
redshift subsample. The typical procedure to convert 
from higher-J CO measurements to Sgas involves cor¬ 
recting the CO luminosity for subthermal excitation to 
arrive at the CO(l-O) equivalent, using aco to infer the 
molecular gas mass, and division by one of the available 
size measurements to get the surface density. However, in 
Section [5.1 1 we found that other measurements of galaxy 
size systematically underestimated the low-J CO size by 
a factor of ~1.3, which corresponds to a 1.7x underesti¬ 
mate in the area and overestimate in Sgas. This error is 
shown as a horizontal arrow in Fig.[^ and is very similar 
to the offset seen between the two nigh-redshift galaxy 
samples. Clearly, care is needed when placing galaxies 
observed in heterogeneous ways on the SF relation. Fur¬ 
ther comprehensive studies of low-J CO emission by the 
VLA or ATCA offer the potential to resolve this issue 
by directly comparing the effective radii of high-redshift 
galaxies at various CO transitions. 
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Figure 7. The Schmidt-Kennicutt star formation r elation. For 
all ga l axies, we use CO-H 2 conversion factors using the |Narayanan| 
|et al.| ||2012|l formula. The high-redshift galaxies are separated into 
those that nave been resolved in a low-J CO transition and those 
for which a higher line was observed (with a size determined from 
other means, generally the rest-frame optical/UV size). The black 
arrow shows the effect of the average correction to Sgas that results 
from the difference in size between the low-J CO emission and the 
rest-frame UV or dust emission, as seen in Fig.[^ The local galaxy 
sample is from 
compi l ed fro m 

l et al.l <| 20Tot: [T _. ■ . ... _. . - ... _ 

(|2U12 i; | Hoclge et a l.J (|2U12 ;|Walter et al .l (|2Ui:^|); Ihu et al.l (|l^U 13l); 
|lvison et al.| l |2Ui3[ |; flacconi et al.| l |^13[ |; [rLodge et ah] ( |2Ul5^ ' 


6. CONCLUSIONS 
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We have presented spatially and spectrally resolved 
images of two gravitationally lensed dusty star-forming 
galaxies at redshifts z = 2.78 and z = 5.66. SPT0346-52 
is among the most intrinsically luminous DSFGs, while 
the luminosity of SPT0538-50 is more typical of the 
DSFG population. Using a visibility-based lens mod¬ 
eling procednre, we have shown that SPT0346-52 has 
complex dynamics, and confirmed the merger hypothe¬ 
sis in SPT0538-50. By comparing with lens models de¬ 
rived from ALMA observations of the dust continuum 
in each galaxy, we find that the difference in magnifica¬ 
tion between the molecular gas and dust varies between 
0 and 50%, mostly due to the larger physical extent of 
the gas compared to the area of active star formation. 
In SPT0538-50, we have shown that the physical extent 
of the CO emission decreases with increasing transition, 
with the CO(3-2) emission being roughly the same size 
as the dust continuum. We have constrained the CO-H 2 
conversion factor via three independent methods, find¬ 
ing values near those expected for highly star-forming 
systems. The three methods agree reasonably well when 
applied to these two objects; further in-depth studies 
may be able to discern systematic differences between 
the various methods. Finally, we have placed these two 
objects on the Schmidt-Kennicutt star formation rela¬ 
tion, finding that they he along the upper envelope of 
vigorously star-forming systems. Part of this offset may 
be explained by the different effective source sizes of the 
CO emission as a function of observed transition, an ef¬ 
fect that should be taken into account as larger samples 
of spatially-resolved high-redshift molecular gas measure¬ 
ments become available. 
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